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We review and compare two different approaches to radiation reaction in classical electrodynamics 
of point charges: a local calculation of the self-force using the charge equation of motion and a 
global calculation consisting in integration of the electromagnetic energy-momentum flux through 
a hypersurface encircling the world-line. Both approaches are complementary and, being combined 
together, give rise to an identity relating the locally and globally computed forces. From this 
identity it follows that the Schott terms in the Abraham force should arise from the bound field 
momentum and can not be introduced by hand as an additional term in the mechanical momentum 
of an accelerated charge. This is in perfect agreement with the results of Dirac and Teitelboim, 
but disagrees with the recent calculation of the bound momentum in the retarded coordinates. We 
perform an independent calculation of the bound electromagnetic momentum and verify explicitly 
that the Schott term is the derivative of the finite part of the bound momentum indeed. The 
failure to obtain the same result using the method of retarded coordinates tentatively lies in an 
inappropriate choice of the integration surface. We also discuss the definition of the delta-function 
on the semi-axis involved in the local calculation of the radiation reaction force and demonstrate 
inconsistency of one recent proposal. 

PACS numbers: 04.20.Jb, 04.50.+h, 46.70.Hg 



I. INTRODUCTION AND OVERVIEW 

Studies of the radiation reaction in classical elec- 
trodynamics initiated by Lorentz and Abraham as 
far as in the 19-th century, remained an area of ac- 
tive research during the whole 20-th century. Al- 
though with the development of quantum electrody- 
namics this problem became somewhat academic, it 
still attracts attention in connection with new appli- 
cations and new ideas in fundamental theory. The 
current understanding of the radiation reaction has 
emerged in the classical works of Dirac 111, Ivanenko 
and Sokolov 21, Rohrlich , Teitelboim Q and some 
others. One of the crucial points of this theory is the 
nature of the non-dissipative term in the reaction force 
known as the Schott term p|. According to Dirac 
and Teitelboim p|, this term is a finite part of the 
derivative of the momentum of the electromagnetic 
field which is bound to the charge. Meanwhile, as far 
as we are aware, in the existing literature there are 
only two explicit calculations of the bound momen- 
tum via integration of the corresponding flux: that 
by Dirac and a more recent calculation by Poisson [(| , 
who used the method of retarded coordinates of New- 
man and Unti 0. The results disagree, namely, ac- 
cording to 0, the bound momentum produces only 
an infinite term which has to be absorbed by the mass 
renormalization. The purpose of this paper is to clar- 
ify the nature of the above discrepancy and to provide 
an independent calculation of the bound momentum. 

In the famous paper published in 1938, Dirac P 
gave a detailed and consistent derivation of the equa- 
tion of motion of a point radiating charge in four- 
dimensional flat space-time. In his original formula- 
tion, the decomposition of the retarded field into the 



sum of the radiation field (the half-difference of the 
retarded and advanced fields) and the self-action field 
(the half-sum of the same quantities) was suggested. 
When such a decomposition is used to describe the 
action of the proper field on the charge itself, the first 
part provides the finite radiation-reaction force (the 
Abraham vector 0) 

/^ = ^(«V + /), (l) 

where is the four-velocity and a M = is the four- 
acceleration. Here the dot denotes the derivative with 
respect to the proper time, a 2 = a^a^ 1 , and the signa- 
ture + ,—,—,— is understood. The second part con- 
tributes only to an infinite renormalization of mass 
and does not affect the motion of the charge in the 
flat space-time. 

The reaction force contains the radiation recoil 
term, the first term in 

/Lit = fe 2 aV, (2) 
and the Schott term 

Schott = (3) 

This second term is a total derivative, so it does not 
correspond to an irreversible loss of momentum by 
the particle, but plays an important role in the mo- 
mentum balance between the radiation and particle 
momentum loss. If / r ^ d = 0, this does not necessarily 
mean that there is no radiation (recall the well-known 
case of radiation from a uniformly accelerated charge 

Hi), 

but if there is no radiation, the Abraham force 
d is zero. Indeed, if one has o M a M = at any time, 
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then it is easy to show that the three-acceleration is 
zero, a = 0, and therefore a M = 0. Thus, no radiation 
reaction force is possible in the absence of radiation. 

The correct interpretation and a constructive 
derivation of the Schott term was given by Dirac him- 
self pj via integration of the Maxwell stress-tensor of 
the charge retarded field over the space-like hyper- 
surface orthogonal to the world-line. The subsequent 
discussion was somehow obscured by the use in this 
paper of the advanced potential. But Havas [Ty] no- 
ticed that actually only the physical retarded field was 
involved in the calculation of the bound electromag- 
netic momentum, and the whole expression for the 
Abraham vector can be obtained using only the re- 
tarded field E2- This became especially transparent 
after a later investigation of the nature of the Schott 
term by Teitelboim |4| (see also the review E2)) where 
it was emphasized that this term originates from the 
bound electromagnetic momentum. However, having 
provided a very clear and comprehensive discussion, 
Teitelboim did not present an explicit calculation of 
the quantities involved, addressing the reader to the 
Dirac's paper for technical details. Meanwhile, the de- 
tails of the integration carried out in are fairly non- 
trivial, and that is why more recently an attempt was 
made in Ref. @ to simplify the derivation using the re- 
tarded coordinates of Newman and Unti ■ This sim- 
plified method seems to be getting popular and it has 
been generalized to arbitrary dimensions [13j. How- 
ever, this modified calculation fails to give the Schott 
term as a part of the bound electromagnetic momen- 
tum. This has led the author of .6] to revive the at- 
tempt to ascribe a mechanical origin to the Schott 
term, which interpretation can in fact be found in the 
earlier literature. 

The crucial point in Dirac's calculation was the 
power series expansion of the retarded field in terms 
of the suitably defined small parameter related to 
the proper time difference between the moments of 
the emission and observation. This is necessary be- 
cause the integrand of the associated integral expres- 
sion contains retarded fields taken at different mo- 
ments of the proper time. Although one might think 
that such expansions can be avoided by using the 
retarded coordinated, a more careful analysis shows 
that this is not so. We present here a straightfor- 
ward calculation which is technically slightly differ- 
ent from that used by Dirac, but is similar concep- 
tually. We have also generalized the calculation to 
arbitrary space-time dimensions (to be presented else- 
where), which sheds new light to the problem of the 
Schott term. Note that multidimensional generaliza- 
tion of the Maxwell theory was discussed by Ivanenko 
and Sokolov [l4|, soon after the work of Dirac, in con- 
nection with the Huygens principle, and recently this 
problem has attracted attention in view of general in- 
terest to space-time models with large extra dimen- 
sions [3 111 Hi E 13). 

Some confusion about the Schott term is also re- 



lated to the well-known 'phenomenological' derivation 
of the Lorentz-Dirac equation, as given in the book 
by Landau and Lifshitz |19|. In this heuristic deriva- 
tion, the first term of the Abraham force is obtained 
by computing the rate of radiation, while the Schott 
term is added by hand from the requirement of or- 
thogonality of the reaction force to the particle four- 
velocity. Formally, this procedure leads to the correct 
equation (though does not answer the question about 
the physical origin of the Schott term), so per se it 
does not contradict to the correct interpretation ac- 
cording to which the Schott term is viewed as a part 
of the bound momentum. But as we show here, if one 
does not relate the Schott term to the electromagnetic 
momentum, the energy-momentum balance equations 
become contradictory (Sect. II). The essential dif- 
ference between the 'phenomenological' derivation of 
the Schott term via the orthogonalization procedure 
and its consistent treatment as the derivative of the fi- 
nite part of the bound electromagnetic momentum be- 
comes especially clear in higher dimensions. It turns 
out that generically the number of possible momen- 
tum 'counterterms' in higher dimensions is larger than 
the number of equations arising from the requirement 
of the orthogonality. As a result, the Schott term(s) 
can not be obtained within the orthogonalization pro- 
cedure in even dimensions higher than six. 

The redefinition of the mechanical momentum of 
a radiating charge would also be wrong conceptually, 
since it would imply that the Maxwell-Lorentz elec- 
trodynamics be not muliplicatively renormalizable: a 
new (finite) mechanical term not present in the initial 
lagrangian would be required. In fact, failure of the 
multiplicative renormalizability is what hap pen s in- 
deed in higher-dimensional electrodynamics |l6t Il7| , 
but an essential property of the four- dimensional the- 
ory is the renormalizability in the sense that no ad- 
ditional counterterms (either infinite or finite) are re- 
quired to make the theory consistent. This reflects the 
renormalizability of quantum electrodynamics in four 
dimensions. 

Our other remark concerns the regularization of 
products of the delta-function and its derivatives with 
the Heaviside function which arise in the quasilocal 
treatment of the radiation reactionproblem. Recently 
this problem was reconsidered in E3 m the context of 
the higher-dimensional generalization of the Maxwell 
electrodynamics, where a new regularization scheme 
was suggested and used to derive the Lorentz-Dirac 
equation in even dimensions higher than four. We 
will show below that the regularization proposed in 
17] fails to reproduce the correct result already in 
the case of four dimensions, so the validity of the 
equations derived in 0, (and of the regulariza- 
tion itself) is questionable. Meanwhile the consistent 
treatment of the delta-function on a half-line via the 
point-splitting was suggested long ago (see e.g. the 
books by Ivanenko and Sokolov Rohrlich [3j and 
Barut [2(|), and we have checked that it works per- 
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fcctly well in any even space-time dimensions. The 
actual check of validity of any regularization of this 
kind, as was emphasized by Dirac, consists in an alter- 
native calculation of the reaction force by integrating 
the rate of the variation of the bound electromagnetic 
momentum. 



II. FIELD DECOMPOSITION VERSUS 
STRESS TENSOR DECOMPOSITION 

Our definitions closely follow those by Rohrlich Q 
and Teitelboim [4|. The retarded potential generated 
by a point charge moving along the world-line x^ = 
z fi (s) depends on the kinematic variables taken at the 
(position dependent) retarded proper time s Te t(x) de- 
fined as the solution to the equation 



0, B? = a? - z»( 8ttit ), 



(4) 



satisfying x° > z°. The second solution to the same 
equation with z° > x° defines the advanced proper 
time s a dvi x )- Introducing the invariant distance 



(5) 



which is equal to the spatial distance |R| = |x— z(s ret )| 
between the points of emission and observation in the 
momentarily co-moving Lorenz frame at the time mo- 
ment x° = z (s re t), one can present the retarded po- 
tential as 



Sret (a:) 



(6) 



It is convenient to introduce the null vector c M = 

W I p, whose scalar product with is equal to unity, 

and also the unit space- like vector — c M — v^. Thus 
we have 



1. 



0, 



1. 



■1. 



(7) 



(Here and below we omit, where unambiguous, brack- 
ets in the four-dimensional scalar products, e.g. vc — 
v^Cfj,). Differentiations with respect to x 11 are per- 
formed using the relations 



duP = v f _ l + Ac p , 



(8) 



where 



A = p = p{ac) - 1. (9) 
Using these formulas we obtain the field tensor 
e(p(ac) - 1) 



F„ 



V[nCv] - -a[^c v ] (10) 



where square brackets denote antisymmctrization 
without factor 1/2, e.g. a^c v ^ = a^c^ — a v c^ and 
all quantities have to be taken at the moment s re t(x) 
of the proper time . Similarly, the advanced poten- 
tial A£ d an d the corresponding field can be written 
in terms of quantities depending on s a dv(x). 

Separation of the radiation from the total field can 
be performed in two different ways. The first proce- 
dure is linear in the field and consists in the splitting 
the retarded potential A^ ct into the radiative part 



^rad ~~ 2 ^ ret advJ 



and the 'self part 



4^ — _ i A^ _|_ A^ \ 
A se\i — 2 ^ ret advi 



(11) 



(12) 



The radiative potential satisfies the homogeneous 
D'Alembert equation and changes the sign under re- 
flection of time, as expected for the radiation irre- 
versibly lost by an accelerated charge. This part of 
the retarded field tends to zero in the static limit. 
The self part is time-symmetric and remains finite in 
the static limit, where it coincides with the Coulomb 
potential. 

The split of the second kind is quadratic in the field 
and uses the energy-momentum tensor 



(13) 



Considering the energy-momentum tensor we will 
always deal with the retarded solution of the 
D'Alembert equation. Substituting here the field ten- 
sor HI 0|l we obtain the sum of two terms 



emit bound 1 



where the first term is proportional to p 



(14) 



47T „,. 
~~^2 emit 



{acf +a 2 



-c^c 



(15) 



while the second one contains higher powers of p 



g 2 bound 



+ 



a^c^ + 2{ac)cfc v - (ac)v^c^ 
? 



(16) 



Here the symbol (pis) denotes a symmetrization with- 
out the factor 1/2. 

The first expression (|15(l is distinguished by the fol- 
lowing properties: 

• its geometric structure is the tensor product of 
two null vectors c M , 

• it is traceless, 

• it falls down as |x| -2 when |x| — > oo, 
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it is independently conserved 



d = 



(17) 



This latter property follows from the differentiation 
rules ©• All these features indicate that T^ it de- 
scribes an outgoing radiation. An independent con- 
servation of this quantity means that the bound part 
is also independently conserved 



u " ± bound — 



(18) 



Conservation of the total four-momentum implies 
that the sum of the mechanical momentum and the 
momentum of the electromagnetic field is constant 
(for simplicity we do not include an external field ): 



nxech 



ds 



d Pcm = q 

ds 



(19) 



Here the mechanical part is proportional to the bare 
mass of the charge 



Pmech = TO 0^, 



(20) 



Now, somewhat unexpectedly, / r ? d 
/self /bound- differing by the Schott term 



^ /emit and 



J rad 

_ ft* 
•'self — Jb( 



f M -4- f M 
•/omit T J Schott' 



bound 



Schotf 



(27) 
(28) 



The identity H26|) is satisfied as expected. Explicit 
calculations verifying these results will be presented 
in what follows. They are fully consistent with the 
results of Dirac 0, Rohrlich Q and Teitelboim Q. 
On the contrary, in it was found that f£ ound = 
/self' while the Eq. J23) still holds. This is obviously 
inconsistent with the energy-momentum conservation 
identity pjjl . 

To avoid any confusion, we note that though both 
sides of the energy-momentum conservation identity 
l|26|) contain divergent terms, the extraction of finite 
terms is fully unambiguous because of their different 
dependence on the kinematical variables. Moreover, 
the parameterization of the divergent terms can be 
made similar in both calculations, so these terms mu- 
tually cancel in (|2"oT) before the regularization is re- 
moved. 



while the field part is given by 



jC = / T^dZ,, 



(21) 



where integration of the electromagnetic stress ten- 
sor is performed over a space-like hypersurface whose 
choice will be discussed in detail later on. It has to 
be emphasized again that in the expression (| 1 31) for 
the stress tensor of the electromagnetic field one has 
to use the physical retarded field. According to the 
second splitting, one can write 



dp" 



mcch 



+ /£ 



./ be 



dPcm rl-i elt 

'da < 
d 

ds 
d 
ds 



ound' 



rpliv jy 
J -cmit UZj i / ' 



J bound UZj '- 



(22) 
(23) 
(24) 



On the other hand, the derivative of the bare me- 
chanical momentum can be expressed using the equa- 
tion of motion of the charge in which the electromag- 
netic field is decomposed into the self part (|T2)l and 
the radiation part (|f f I) 



dp? 



mcch 



ds 



_ ft , ft* 

J self J rad ' 



(25) 



Clearly, the following energy-momentum conservation 
identity should hold in view of l|2"2")) and (J2EJ): 



f + f M = f M 4- f 

•/self 1 -/rad -/bound 1 -/emit' 



(26) 



III. WORLD-LINE CALCULATION: POINT 
SPLITTING 

The retarded and advanced potential taken on the 
world- line = z M (s) of a charge can conveniently be 
written in terms of Green's functions 0] 



G seU (Z) = 6(Z 2 ), 
G T3A {Z) = i^S(Z 2 ), 



(29) 
(30) 



where Z» = Z^(s,s') = z^{s) - z^(s'). Substitution 
of the electromagnetic field of the charge on its world- 
line leads to the following integrals 

r(s)=2e 2 J Z^(s,s')v^(s')v v (s)-^G(Z)ds\ 

(31) 

for both and f{*,. Due to the presence of delta- 
functions in G sc if and G ra( j, one is tempted to expand 
the integrands in a = s — s'. Taking into account that 
Z 2 = a 2 + 0(<t 4 ), one can write 

G seli (Z) = <5(a 2 ) + 0(cT 4 ), (32) 
G rad (Z) = ^(J(a 2 ) + 0(cr 4 )). (33) 



Expanding the rest of the integrands in tr, one encoun- 
ters the following integrals: 



r°° d 

J-oo da z 



(34) 



for the self-force, and 
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for the radiation reaction force, with I > 2. 

Both these integrals are ill-defined. Passing to the 



consists in the 'point-splitting' 



variable x = 
of the type 



a they can be transformed to integrals 



(f>{x)S [k) (x)dx, 



(36) 



where S ^ (x) denotes the k-s derivative of the delta- 
function, which function has the support at the 
boundary point of the integration domain. Comput- 
ing such an integral is equivalent to taking the prod- 
uct of the delta-function and the Heaviside function, 
or, equivalently, defining 5(x) on the semi-axis. For 
this a suitable regularization is needed. Before dis- 
cussing this point, we note that, with any regulariza- 
tion, the integrals l|34|) . (|35[1 should vanish for I > 3 
by power counting. Moreover, all Ai vanish for odd 
I, and all Bi vanish for even I from parity consider- 
ations. By power counting one can also show that 
terms in (13211 , (|33|l proportional to er 4 actually give no 
contribution, while the relevant terms that must be 
retained in the expansion of the integrand are given 

by 

2X^(s, s'KVK(s) = a»a 2 -\{^ + v^a 2 )a 3 . 

(37) 

Now we discuss the meaning of the delta-function 
on the semi-axis. Recently an attempt was made 17] 
to develop a general theory for such objects using pri- 
mary definitions of the theory of distributions. The 
proposed regularization for the first derivative of the 
delta-function defined on the semi-axis (see Eq. (23) 
in 0) reads: 



5'(x) = lim 



a^+a da a ' 
Using this regularization one finds: 



x > 0. 



(38) 



6(a 2 ) = lim S(a 2 -e 2 ) 



lim 



S(a - s) + 6 (a + e) 
2e 

with a prescription that the limit should be taken only 
after evaluating all the integrals. With this, omitting 
the symbol of the limit, we obtain 



A-2 

B, 



1 

2? 
-1, 



so one finds 



J self 



2s 
2J 2 
3 



After the mass renormalization, 

too - A2 = m, 
we get the Lorentz-Dirac equation 
2e 2 

ma^ = (v^a 2 

3 V 



(44) 
(45) 

(46) 
(47) 

(48) 
(49) 



One has to realize, however, that there is no purely 
mathematical proof of the correctness of the point 
splitting precedure. So actually to check the validity 
of this or any other regularization involved in the local 
calculation of the radiation reaction force one has to 
perform an alternative calculation via the integration 
of the electromagnetic stress-tensor. 



IV. INTEGRATION OF 
ELECTROMAGNETIC MOMENTUM 



A. General setting 



A-? = lim 



4-y/a' 



5, = 1. 



(39) 
(40) 



Substituting this into the above formulas we obtain 



A* _ 



J self 



J rad 



lim — - — " M 



2e 2 



{v^a 2 + d"). 



(41) 
(42) 



The first quantity can be absorbed by the renormaliza- 
tion of mass, while the second gives the correct struc- 
ture for the Abraham force but with the wrong sign. 
Thus the proposal f° r the delta-function and its 
derivatives on the semi-axis fails to give the correct 
result. 

Meanwhile, a satisfactory method to deal with such 
integrals was suggested long ago (see 0, m H3)i ^ 



An alternative derivation of the radiation reac- 
tion force consists in the integration of the energy- 
momentum tensor of the electromagnetic field of a 
charge over an appropriate hypersurface in space- 
time. The correct choice of this hypersurface is es- 
sential for the calculation, so we will discuss it in de- 
tail. Assuming the split (|14|l of the stress-tensor into 
the emitted and bound parts, one can consider the 
emitted and bound momenta separately. Teitelboim 
||j defined the corresponding integration surfaces dif- 
ferently in both cases taking into account the specific 
properties of these quantities. Here we will show that 
both integrals can be transformed to those over an in- 
finitely thin world tube around the particle world-line. 

We would like to calculate the four-momentum car- 
ried by the electromagnetic field of the charge for a 
given moment of the proper time s on the particle 
world- line z^(s). To do this one has to choose a space- 
like hypersurface S(s) intersecting the world-line at 
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FIG. 1: Integration of the bound electromagnetic momen- 
tum. Here £(si) is the space-like hyperplane transverse to 
the world- line z fl (s) intersecting it at the proper time si 
(similarly E(s2)). The hypersurfaces S e and Sr are small 
and large tubes around the world- line formed by sequences 
of the 2-spheres dY e (s) and 8Yr(s) for s 6 [si, sz]. The 
domain Y(s2) C £(s2) (similarly Y(si)) is the 3-annulus 
between 8Yr(s2) and dY e (s2). 



the point z M (s) and to integrate the electromagnetic 
energy-momentum flux as follows 



rem 



00 



(50) 



E( S ) 



The simplest and the most practical choice for S(s) is 
that of the hypersurface orthogonal to the world-line 



(51) 



The integral (|50|) is divergent near the world line. 
The rate of the divergence can be controlled introduc- 
ing the small length parameter e, the radius of the 
2-sphere dY e (s) (Fig. 1), defined by the intersection 
of the hyperplane (|51() with the hyperboloid 



z{s)f 



(52) 



We also introduce the sphere dY R {s) of a large radius 
R defined by the intersection of S(s) with the hyper- 
boloid 



-R 2 



(53) 



The electromagnetic momentum can then be obtained 
by taking the limit e — » 0, R — > oo of the integral 
over the domain Y(s) C S(s) between the boundaries 
dY e (s) and aifl(s). 

Let us evaluate the variation of this quantity be- 
tween the moments Si and S2 of the proper time on 
the world-line of the charge 



At/ 1 



(54) 



Y(s 2 ) 



For the bound momentum it is convenient to con- 
sider the tubes S E and Sr formed as sequences of the 
spheres dY e (s) and 8Yr(s) on the interval s € [si, sz] 
and to transform this quantity to 



Mo 



-'bound" " 



-'bound" " 



(55) 



in view of the conservation equation for T^ und JT^J. 
Here normal vectors in dS v are directed outwards 
with respect to the world-line. The contribution from 
the infinitely distant surface Sr vanishes if one as- 
sumes that the charge acceleration is zero in the limit 
s — > — oo 4j. This assertion is somewhat non-trivial, 
since, in spite of the fact that the stress tensor (|16|l 
decays as R~ 3 , the corresponding flux does not van- 
ish a priori, because the surface element contains a 
term (proportional to the acceleration) which asymp- 
totically grows as R 3 (see the Eq. EH1 below). As a 
consequence, the surviving term will be proportional 
to the acceleration taken at the moment s ro t of the 
proper time, where s rc t — ► — oo in the limit R — > oo. 
Finally we are left with the integral over the inner 
boundary only 



Abound 



boun 



^dSy. 



(56) 



For integration of the emitted momentum it is con- 
venient to take the light cone boundary C(s') instead 
of Sr as shown in Fig. 2. The actual change of the 
emitted momentum in the whole three-space corre- 
sponds to the limits s' — > — oo,e — > 0. Since the 
normal to the light cone lies on it, the flux of the 
energy-momentum tensor through the null boundary 
C(s', Si, S2) is zero for any s', therefore 



Z(s',s 2 ) 



Z(s',s 1 ) 



emit 

S s C(s',S2,Sl) 



Tcmit C vdS 



rpfiU ,q 



(57) 



Se 
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dZ(s',s 2 ) 




FIG. 2: Integration of the emitted momentum. Here C(s') 
is the future light cone of some point s' on the world- 
line, and C(s',si,S2) is its part between the hypersur- 
faces and £(52). The domain Z(s',S2) (similarly 

Z(s',si)) is the annulus between the intersection of the 
light cone with E(s2) (the outer boundary) and the small 
sphere 9Y e (S2) (the inner boundary). 



where the limits e — ► 0, s' — > —00 have to be taken. So 
both the emitted and bound momenta can be reduced 
to integrals over the small tube around the world-line. 

Now we have to find an integration measure on the 
small tube. 



B. Induced metric on S e 

Consider the foliation of the space-time region 
shown in Fig. 1 by the hypersurfaces E(s) parameter- 
ized by the spherical coordinates r, 9\ = 9, 9 2 = tp. In- 
troducing the unit space-like vector n^(s,9i), n^n^ = 
— 1, transverse to we can write the following coor- 
dinate transformation from x^ to the set s,r,0i\ 



x» = z' M (s)+rn fl (s,e i ), 

^( s )/(Mi) = o, 



(58) 
(59) 
(60) 



Since the 4-acceleration vector a M is orthogonal to the 
4-velocity, it lies in the hyperplane £(s) and we can 
further specify the angular coordinates choosing the 
polar axis along the three-acceleration a. Our conven- 
tion about the angle variables is such that the four- 
vector in the Lorentz frame momentarily comoving 
with the charge at the proper time moment s, CF(s), 
is given by 



Then a^n^ = —acos9, where a = |a|. Differentiating 
Ij58(l with respect to the new coordinates we obtain: 



dx^ 
~d9~ 



dn»{s,9i) 



d9i ' 



(62) 
(63) 
(64) 



Let us calculate the derivative dn^{s,9i)/ds in 
CF(s). To find the variation n^(s + ds) — n%, we 
consider another Lorentz frame, CF(s + ds), which is 
comoving with the charge at s + ds. These frames 
are related by the Lorentz boost with the velocity 
|v(s + ds) — v(s)| = |a(s)ds| = ads. Performing the 
Lorentz transformations and taking into account that 
cosf? is the same in both frames one finds n^{s + ds) 
in CF(s): 

tt/' (s + ds) — (acos9ds, sin 9 cos ip, sin 9 sin ip, cos6>) . 

Hence in this frame 

rt l (s + ds) - n»(s) 



ds 



(a cos 9,0, 0,0), 



which is proportional to the four-velocity in the same 
frame v IJ '=(l, 0,0,0). Thus in an arbitrary frame one 
has 



— — = acos#v M = -(an)/. 

OS 

Denoting the angle derivatives as 
e£=rdni>(s,e i )/d6i, 



(65) 



wc find the relations 



^ M eJ = r 



dn 2 
'W9~ 3 
d(vn) 
d9i 



rn »d9- = °> 



(66) 



n^v H 



0. 



showing that the four vectors {n^, e^ 1 } represent 
the space-time vierbein, with being the tangent 
vectors to the 2-sphere. The induced metric on the 
sphere is gij — •q^e^e 1 ^, and for the four-dimensional 
induced metric we obtain 



9^ 




and therefore 



dct g = [1 — r(an)] 2 detgij. 



(67) 



(68) 



Hence the area element on the hypersurface r =const 
S s will be given by 



(0, sin 6* cos ip, sin# sin ip, cos 9). (61) 



d 3 a = r 2 [l - r(an)]dsd£l, 



(69) 
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where dfl = sm.6d9<kp. Setting r = e one finds for the 
induced metric on S F 



d 3 (j = e 2 [l - e(an)]dsdfl, 



(70) 



From Ijtitifl one can see that n M is a unit vector normal 
to the tube, so finally 



dS^ = e 2 [\ — e(an)]n IJ ,dsdVL. 



(71) 



Substituting this into I|57I56(I we find for the deriva- 
tive of the electromagnetic momentum the following 
expression 



ds 



J e 2 [l - e(an)]T£fa,<Kl, (72) 



valid for both the emitted and bound parts, where 
the limit e — » is understood. It has to be realized 
that the dependence of T££ on e is somewhat non- 
trivial. In fact, the energy-momentum tensor depends 
on the space-time point through the quantity p, de- 
pending directly on x M , and also through the retarded 
proper time s rc t- At the same time, we need to express 
the resulting quantity as a function of the proper time 
s corresponding to the intersection of the world-line 
with the space-like hypersurface. Carefully keeping 
track of all this we expand the stress-tensor in terms 
of e as follows 



Tfiv (p : S r 



]T e fe e^( s ,^,^, a ^,^). 

k=-4 



where all quantities in the right hand sides are taken 
at the proper time moment s. Omitting also the e- 
term in the integration measure in (|57J) and assuming 
Si = — oo, S2 = s we obtain 



Pemit( S ) 



4tt 



ds' (a 2 + (an) 2 ) («" + n")dfi. 



(76) 



After integration over the angles we arrive at 

s 

— oo 

so the 'emitted' contribution to the radiative force is 

dn^ 2e 2 
/lit - "%* = 4"aV. (78) 



D. Bound momentum 

In this case calculations are substantially more in- 
volved. All quantities in i|I6[) depend on the retarded 
time, and to facilitate their expansion in e-series it 
is useful to express T^° und through the null vector 
Rt* = c f *p: 

4tt t ^ = a^R v) (2(aR)-l)R»R v 

g 2 1 bound p 4 p 6 

(1 - (aR)) v^R^ rf v 

+ ~p~ 5 V' 

The expansion of i?' 1 reads 



(79) 



The leading divergent terms here are proportional to 

e~ 2 for the emitted part (QH) and to e" 4 for the bound rm = x n _ z ^( Srct ) = x » - z^(s) + z M (s) - ^(s rct ) 



part JT5J. Substituting this expansion into ((72)) we 
have to perform integration over the angles and then 
to take the limit e — *• 0. The integration is easily done 
using the formula 



47T 

n^rivdVL = — Apt,, A M „ = v^v v - g^, (73) 



while the integration of products of an odd number of 
n^ gives zero. 



= en* + xPo - ^a 2 + id^cr 3 + C(a 4 ), (80) 

where a = s — s rc t > and all the vectors in the last 
line are taken at s. This is an expansion in powers of a, 
but we need an expansion in powers of e. The relation 
between the two can be found from the condition R 2 = 
0. Assuming an expansion of a in terms of e 



<7 = ^& fc £ fe , 



(81) 



C. Emitted momentum 

Considering first the emitted momentum case, we 
see from Eq. (|15|l that the area factor e 2 compensates 
the denominator p 2 (e) ~ e 2 in the expression for T^ it , 
so in the limit e — > it is sufficient to take only the 
leading terms in the numerator 



and substituting it into the equation R 2 — 0, one finds 
order by order the coefficients: &o = 0, b\ = 1, b 2 = 
an/2 and so on. Thus, up to the third order terms, 
which is sufficient for our purposes, we obtain: 

an s 

cr = e -\ e 2 + (9(an) 2 + a 2 - 4an) — . (82) 

2 v 1 24 

Substituting this into Eq. H80[l we find: 



I S rct 



(ac) 



a 

n h + v 

2 



(an) , 



(74) 
(75) 



RP= (n fl +v fl )s + ((an)v tl -a^)—+ (83) 

£ 3 

+ [(9(an) 2 + a 2 - 4dn) - 12(an)a M + 4d M ] — . 
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Similar expansions for the velocity and the accelera- 
tion taken at the moment s ret are 



v^\ StBt =v li - a^e + {a» - (071)0") —- (84) 



[(9(an) 2 + a 2 - 4dn) a^ - 12(an)d^ + 4a>*] 



24' 



a"| Srct = - d"e + (a" - (an)/) —. (85) 

The invariant distance parameter p — u /J (s re t)-R' 1 is 
given by the product of the two expansions: 

2 3 

p = e- an%^ + (8(dn) - 3a 2 - 3(an) 2 ) |-. (86) 

For convenience we give also the expansion of the most 
singular term in (|79|l up to the relevant order: 

1 1 r 

— = — 1 + 3 {an)e+ 
p 6 e 6 

+ (6(an) 2 + 3a 2 /4 - 2dn) e 2 } , (87) 

and the expansion of the scalar product (aR): 

aR = an + (a 2 /2 — an) e + [4da/3 + a 2 (an)+ 
+dv + an - (an) (an)] e 2 /2. (88) 

Substituting all these expansions into (|56|l we obtain 



AK ound = ^ / ds \^ + ^ + [((an) 2 + a 2 /3) ^+ 



+ ((an) 2 + a 2 /2) n" - 2a^/3 + 3(an)a fi /4] } dfi. 

(89) 

Using the above rules of the integration over the an- 
gles l|73f) . one can see that the leading divergent term 
proportional to 1/e 2 vanishes and the result reads 



Ap 



bound 



/ ,/. s f JL m - ^d" 



(90) 



Setting si — — oo, S2 = s, we obtain 

^-= e2 /l ds '(i aM -H- (9i) 

Therefore the bound part of the self-force is 



J bound 



^ou„d e 2 a^ 2e 2 A 



2e ^ 



Here the first divergent term has to be absorbed by 
the renormalization of mass, while the second is the 
Schott term. Comparing this with 1461) we confirm 
the relation igHJ and the identity The Schott 

term therefore is the finite part of the derivative of the 
bound electromagnetic momentum of a charge. Note 
that a priori the parameter of regularization e (the 



radius of the small tube) is not related to the split- 
ting parameter of the delta- function in the local force 
calculation. But actually they give the same form for 
the divergent term, for which reason we denoted them 
similarly. With this convention, the divergent terms 
in the momentum conservation identity (|26|l mutually 
cancel. 



V. CONCLUSION 

The goal of this paper was to clarify the recent dis- 
crepancy in calculations of the bound momentum of 
the radiating charge and to present an independent 
explicit calculation revealing the nature of the Schott 
term in the Lorcntz-Dirac equation. Wc would like to 
emphasize the role of the global momentum conserva- 
tion in understanding the origin of the Schott term. 
From the momentum conservation one can derive an 
identity relating the contributions to the reaction force 
in two alternative calculations: a local computation 
on the world-line and the global integration of the 
Maxwell stress-tensor. This identity demands that the 
Schott term should arise from the bound field momen- 
tum, and this is confirmed by an explicit calculation. 
Although this problem was discussed by a number of 
authors in the past, only few of them presented the 
details of calculation of the bound momentum, most 
notably Dirac 0, with whom we perfectly agree. We 
disagree, however, with more recent results of Poisson 
@, who used a different method, based on retarded 
coordinated, in order to avoid somewhat lengthy ex- 
pansions used in the Dirac's calculation. According 
to |6J, the bound momentum reduces entirely to the 
mass-renormalization term. As we have shown here, 
this result is incompatible with the momentum con- 
servation identity and thus can not be correct. More 
detailed analysis shows that the integration surface in 
0] does not correspond to the definition of the bound 
momentum associated with the given moment of the 
proper time of a charge (see the Appendix below). 
Therefore we confirm the result of Dirac and show 
that no alternative interpretation of the Schott term 
is needed (and not possible, moreover). Technically, 
our derivation is slightly different from that of Dirac in 
performing all the necessary expansions in a uniform 
way. Conceptually, we fully agree with Teitelboim Q , 
presenting in addition full details of the global deriva- 
tion of the Schott term. 

Notice that in the literature one often encounters 
a confusing notation related to two alternative ways 
to select radiation. Namely, one has to distinguish 
the 'radiative' field as the half-difference of the re- 
tarded and advanced fields and the 'radiative' com- 
ponent of the Maxwell stress-tensor, which is defined 
through the retarded field. Here we suggest a differ- 
ent wording for the radiative part of the stress tensor, 
which we call 'emitted'. This distinction is necessary 
indeed since the corresponding local forces differ by 
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the Schott term. 

Physically, the Schott term describes the reversible 
variation of the Coulomb field bound momentum dur- 
ing an accelerated motion of the charge. Variation 
of the mechanical momentum of the charge consists 
of this reversible part (which may be both positive 
or negative) and an irreversible loss due to radiation. 
In other words, the momentum radiated away can be 
borrowed both from the mechanical momentum and 
from the bound Coulomb momentum, and this ex- 
plains how in the threshold case of the uniform ac- 
celeration the mechanical momentum of the radiating 
charge may remain constant. One is not allowed to 
redefine the mechanical momentum of the charge by 
adding to it the Schott term without facing contradic- 
tion with the global energy-momentum conservation 
identity. 
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A. Appendix. Integration of flux in retarded 
coordinates 

Here for comparison we reproduce the formulation 
of the problem in the retarded coordinates (re- 
cently generalized to arbitrary dimensions in |l3|). 
The main difference with our procedure is the defi- 
nition of the integration hypersurface for the Maxwell 
stress-tensor in terms of the retarded coordinates due 
to Newman and Unti (s = s rc t, p, 0i). These are in- 
troduced as follows. Selecting an arbitrary point z(s) 
on the world-line one constructs its future light cone 
(Fig. 4) 

C{s) = {x\(x-z(s)) 2 = 0,x°^z }, (93) 

ascribing to it the unique value of the coordinate s, 
and defines on it the coordinate p as an affine length 
parameter on a null geodesic specified by the angle 
coordinates 6 1 .;. The orbit of the constant p(s, x) on the 
cone <|93[) forms a two-dimensional manifold E(s,p) 
which is an intersection of the hyperplanc 

H(s,p) = {x\(x - z(s)Yv^(s) =p = const} (94) 

and the cone i|9"5|) . The open tube surrounding the 
world- line is then defined as the sequence £(it, p) of 
hypersurfaces p ~ const for u G (— oo,s]. The area 
element on this tube is equal to 




FIG. 3: Integration tube in the retarded coordinates. Here 
C(s) is a future light cone with an apex at z(s), H(s,p) is 
a hyperplane transverse to v(s) which intersects the world- 
line at a point s. The two-dimensional surface S(s,p) is 
an intersection E(s,p) = C(s) f) H (s, p), S p is a lateral 
hypersurface of constant p, formed by all E(s',p), s' = 

— OO...S. 

where cy, = (x — z(it)) M / p. The factor p here is the 
same as in the expressions ()15|l . (I16JI . so it partially 
compensates the denominator of T^ und and fully com- 
pensates that of TgV t , which looks as simplification. 
The procedure proposed in [6j consists in the inte- 
gration of the flux over the tube for fixed p with the 
subsequent limit p — ► 0. 

However, this procedure has a serious drawback 
when applied to the bound part of the electromag- 
netic momentum because of the singular nature of the 
integrand at p = 0. In fact, this means that we have 
to consider the sequence of tubes of variable radius p. 
But with variable p the integral 

s 

J du J T^p 2 (v^ + {p(ac) - l) C/1 )dfi (96) 

— oo 

does not give the field momentum associated with any 
given moment of the proper time. Indeed, when p is 
changing, the spheres E(s, p) move across the light 
cone and do not lie on a definite space-like hypersur- 
face. For any fixed finite p the integral (|96|l is per- 
formed over the hyperplane intersecting the world-line 
at the point of proper time 



dS^ = p 2 {v^{u) + X(u, x)c M (u, x))dudQ, (95) 



s = s + t(s, p) > s, 



(97) 
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as shown on Fig. 3. Only in the limit s — > s one ac- 
tually integrates over the hyperplane intersecting the 
world-line at s, but performing the limit as p — > one 
passes through the sequence of different space-like hy- 
perplanes. Contrary to this, our procedure consists in 
fixing the unique space-like hyperplane and carefully 
expanding the integrand in powers of the radius of the 
tube. This allows to pass consistently to the limit of 
the vanishing tube radius. 

Therefore it is not surprising that a calculation 
along the lines of [|| gives for the bound momentum 



only the leading divergent term (see the Eq. (8.3) of 

B) 

Abound e2 a M 

but fails to produce the second finite term in 192|) . For 
the emitted momentum the procedure of [6j works, 
since the answer in this case is given entirely by the 
leading term. 
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